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Abstract 

There is a new set of variables, which from our analysis appear 
naturally adapted to a Hilbert space description on the reduced phase 
space for vacuum GR. In this paper we present the Hilbert space, 
which features a basis of coherent-like states labelled by the algebraic 
classification of the corresponding spacetime. These wavefunctions, 
which satisfy the semiclassical-quantum correspondence, correspond 
to the solution of the quantum Hamiltonian constraint on the space 
of gauge- invariant, diffeomorphism invariant states in these variables, 
and are free of field-theoretical singularities. 



1 



1 Introduction 



Metric general relativity is a totally constrained system whose action can 
be written as a canonical one-form minus a linear combination of first class 
constraints smeared by auxilliary fields 



Seh[9] = JdtJ^ d^x{7r'^hij -NH- N^Hi) , (1) 

where the phase space variables are the 3-metric on a three dimensional 
spatial hypersurface S and its conjugate momentum, given by (/i^j, tt*-^). In 
equation (1) is the diffeomorphism constraint, given by 



= D,Tr> ~ 0, (2) 

where N"^ is the shift vector and Di is the three dimensional covariant deriva- 
tive with respect to S. is the Hamiltonian constraint, given by 



H = Gijkin'^Tr''^ - Vh R 0. (3) 
N is the lapse function and Gij^i is the metric on superspace, given by 



Gijki = ^ii^jk ~ hijhki). (4) 

Equations (2) and (3) signify the presence of 8 unphysical degrees of freedom, 
which if eliminated would reduce the phase space variables (/ijj, tt*-^) to the 4 
physical degrees of freedom needed to describe GR.^ However, the constraint 
(3) remains unsolved in the full theory due to its nonpolynomial structure 
in the basic variables. 

The constraints can be rewritten in terms of new variables at- 
tributed to Abhay Ashtekar, by canonical transformation from the ADM 
triad description of gravity [1],[2],[3]. In the Ashtekar variables the con- 
straints in smeared form are given by 



H, 



'^The reduced phase space description would entail factoring out the gauge orbits gen- 
erated by the constraints, on the constraint surface. 
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for the diffeomorphism constraint and 

H\N] = ^ d3xiV(^e,,fc6'^^'=5^a^5,^ + e^jke^^at^iB'^) (6) 

for the Hamiltonian constraint with cosmological constant A, where N_ = 
N/ \/ detS is the lapse density function and the shift vector N'' takes on the 
same meaning as in (1).^ The Gauss' law constraint is given by 

Ga[A-o] = J^'^A^oDiK (7) 

which signifies invariance under S'C/(2)_ rotations. Aq is the SU{2) rotation 
angle with covariant derivative Di = {Di)ab = Sabdi + fabcA-l and structure 
constants fabc- 

Equations (5) and (6) can be regarded as the direct analogues of (2) 
and (3) , with the latter being simplified due to the polynomial nature of the 
constraints. However, the presence of an additional constraint (7) entails the 
existence of additional degrees of freedom which are unphysical. This may 
in some ways be seen as a tradeoff for the simplifications in the constraints 
if it is more difficult to eliminate these degrees of freedom in congruity with 
the remaining constraints.^ Therefore, our first order of business will be 
to eliminate the superfluous degrees of freedom in general relativity, which 
entails the introduction of a new set of variables. 



2 Transformation and solution to the constraints 
in Soo/CDJ variables 

We would like to preserve the simplicity of the original Ashtekar variables, 
while at the same time reducing general relativity to its basic physical de- 
grees of freedom. First make the substitution 

= "^aeBi (8) 

^The auxilliary variables N^{x) and A" (a;) can be thought of as the parameters for 
spatial diffeomorphisms and SU{2)- gauge transformations within E, and the auxilliary 
variable N_ can be thought of as the parameter for deformations normal to E. 

'''In loop quantum gravity, the spin network states solve the Gauss' law constraint by 
construction. However, the Hamiltonian constraint remains unsolved in the full theory to 
the present author's knowledge, which leaves open the issue of dynamics. This notwith- 
standing, the loop approach has led to significant physical insights at the kinematical level 
of the phase space. 
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in order to eliminate tlic dcnsitizcd triad 5^ as a basic variable. Equation (8), 
known as the CDJ Ansatz, for the purposes of this paper will be regarded 
as a mathematical identity defining vacuum GR for nondegenerate and 
nondegenerate *ae- The CDJ matrix can further be parametrized by its 
symmetric and its antisymmetric parts 



*ae = CaedV'd + Ke = ^aedlpd + XfOfaOfe, (9) 

where Oae = {e^ '^)ae is a complex orthogonal matrix parametrized by three 
complex angles 6 = {6^ ,9"^ ,6^) and T = (Ti, T2, Ts) are the generators of the 
■so(3, C) algebra in the adjoint representation. Hence, the symmetric part 
Xae of the CDJ matrix with eigenvalues A/ = (Ai, A2, A2) can be seen as the 
result of a Lorentz transformation 



\ae = {e'-')afX}{e-'')fe (10) 

from the frame of eigenvalues, where ^ = 0, into an arbitrary frame parametrized 
by ^ 7^ 0. Note that Im[ff\ and Re[^ correspond to rotations and boosts, 
respectively, in 50(3,(7) language. 

Let us now revisit the constraints in the language of the CDJ matrix 
with a view toward extracting the essential degrees of freedom at the classical 
level. Prom (5), the smeared diffeomorphism constraint is given by 



Hi[N'] = SxeijkN'BiBl^ae- (H) 

For detB 7^ 0, the inverse of the Ashtekar magnetic field exists and 

(11) can be written in unsmeared form as 



{deiB){B-')'lijd = Vx G E (12) 

where we have used (9). The smeared Hamiltonian constraint (6) can be 
written in the form 



H[N] = j d^xN^j^^{yar^ + Adet^') (13) 

where Var'^ = (tr^')^ — tr^'^. While (13) is nonpolynomial in these vari- 
ables, we will see that it is convenient to extract the polynomial part when 
quantizing the theory. Since we are restricting to nondegenerate configura- 
tions detS 7^ and det* 7^ 0, we can then focus on the part 
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Var-^ + Adet* = Vx G S 



(14) 



in solving the constraint. We now compute the ingredients of the Hamil- 
tonian constraint using the parametrization (9). For the determinant we 
obtain 

= det(Aae) + det(eaedAd) + ^eabc^efgl^cgdKehf^d + %fd(^cgd'Ke'4^d'^d>^-{'^^) 

Using the fact that the determinant of an antisymmetric matrix of odd 
rank vanishes in combination with the annihilation of antisymmetric on 
symmetric indices, we end up with 

det^' = detA + -(ea6ce/6d)(e 

efg^ed'g ae 

= detA + ^Xaei^af^cd " 5ad^cf){^ec5 fd' - 5ed'<5/c)V-'(i'0(i' 

= detA + XaeVa'ipe (16) 

where we have made use of epsilon symbol identities. Likewise, we compute 
the variance 

Var^ = {ilXf - (A^e + ^aed'4^d){Ka " ^aed'i^d') = VarX - 2dae1paA- (17) 

The polynomial part of the Hamiltonian constraint (13), which we will from 
now on refer to as the Hamiltonian constraint, then is given by 

H = VarX + AdetA + (AAae - 2Sae)ipai'e = Vx G E. (18) 
Upon substitution of the parametrization (10), then (18) reduces to 

H=^ + ^ + ^ + ^ + {XfOfaOfe - 2Sae)^aA = 0- (19) 
Ai A2 A3 Z 

Note that the complex orthogonal matrix Oae has cancelled out from the 
first two terms of the Hamiltonian constraint, which depend only on the 
invariants of Xae- The last term, quadratic in ipd, depends explicitly on Oae 
which in turn depends on the SO{3, C) frame through 6. 

We would rather like to interpret the Hamiltonian constraint as being 
independent of the SO{3,C) frame and consider the angles 6 as not being 
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independent physical degrees of freedom. The most direct way to do this is 
to recognize that on the space of solutions to the diffcomorphism constraint 
(12), we must have that ^/^^ = 0. For ipd = the last two terms of (19) 
vanish, which enables one to write A3 explicitly as a simple function of Ai 
and A2 

A3 = -f^ ^ ). (20) 

V|AiA2 + Ai + A2^ ^ ^ 

We will ultimately regard the two eigenvalues Ai and A2 as corresponding to 
two physical degrees of freedom for general relativity in the new variables. 

One can see that (12) and (20) constitute a general solution, for non- 
degenerate S^, for the constraints of GR corresponding to invariance under 
spatial and timclike diffeomorphisms. Equations (12) and (20) can be seen 
as the direct analogues of the solutions to (2) and (3), which would suffice to 
construct a solution if not for the Gauss' law constraint. We now argue that 
the Gauss' law constraint is superfluous as follows. Under the CDJ Ansatz, 
(7) reduces to 

Ga[A^]= J^d'xw.i^ae}, (21) 

where we have defined the 'twisted' vector field Wg by 

^ei^ae} = M'^ae} + {fabfSge + /ebg^a/) A^^e^/p (22) 

with Va = Bl^di. The unsmeared form of the Gauss' law constraint is given 
by We{^'ae} = 0, which in the parametrization (9) yields 



WWelV'd} + We{XfOfaOfe} = 0. (23) 

On the space of solutions to the diffeomorphism constraint tp^ = 0, equation 
(23) reduces to 

We{A/(e^^-^)^„(e^^-^);e} = 0. (24) 

Equation (24) is a set of three simultaneous first order partial differen- 
tial equations in six unknowns A^- and 6^ for / = 1,2,3. For each A = 
(Ai, A2, A3), one should in principle be able to find 9 = 9[X; Bf] /9], where [3 
are the boundary data for 6*.^ The point is that while the Gauss' law con- 
straint constitutes a genuine reduction in degrees of freedom of the CDJ 

''The solution to (24) is the topic of a separate paper, which we do not display here. 
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matrix ^ae, it is superfluous when one restricts oneself to the space of 
eigenvalues Xj, since these eigenvalues subject to (20) are freely specifi- 
able. Therefore equation (20) does in fact constitute a reduction of GR to 
two unconstrained degrees of freedom, the physical interpretation of which 
according to (24) fixes a unique SO{3,C) frame 0[X;Bf;/3] by a nonlin- 
ear Lorentz transformation. The implication is that on the gauge-invariant, 
diffeomorphism-invariant subspace of GR, one still has the freedom to choose 
three free functions Ai(x), \2{x) and A3(x). We will exploit this observation 
to describe the dynamics of GR as a constrained system containing only one 
constraint, namely the Hamiltonian constraint, and devote the remainder of 
the paper to performing a consistency check on this dynamics. 



3 Classical equations of motion 

Usually in a theory of first class constrained systems one starts from a phase 
space Lo posessing a cotangent bundle structure and then obtains the reduced 
phase space Ured upon solving the constraints and factoring out the gauge 
orbits. Our approach might seem counterintuitive in that we perform the 
steps in reverse, namely by first imposing the cotangent bundle structure on 
the unconstrained degrees of freedom, while retaining the option to enlarge 
the phase space the constrained version of full general relativity by adding 
in the unphysical degrees of freedom by hand. 

Associate to the eigenvalues of the CD J matrix '^ae, a set of momen- 
tum space dynamical variables ^/ = (*i, ^'2, ^3) and define a phase space 
uj = (X-^,*/), where = {X^,X^,X^) are the set variables on the con- 
figuration space r canonically conjugate to The effect is to impose the 
symplectic structure 

n = (hG)-^ J (fxd-^f{x) A 6Xf {x) (25) 

on the theory at the level prior to imposition of the Hamiltonian constraint, 
with the kinematic constraints having already being taken into account. The 
symplectic two form Q implies the following elementary Poisson bracket 

{Xf{x),^f{y)} = zG<5|<5(3)(x,y) (26) 

for the basic variables. We now proceed to the dynamics of the theory, by 
considering its Hamiltonian evolution at the level of the gauge-invariant, 
diffeomorphism invariant subspace of full GR. Starting from a theory with 
Hamiltonian density 
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H[N] = d^.NideiB)n^,^,^,)y^ + i_ + i_ + i_) (27) 

we will need to obtain the Hamilton's equations of motion, which entails 
the calculation of variational derivatives of (27) with respect to the phase 
space variables {Xf,"^f). The variational derivative with respect to the 
momentum is given by 



6H[N] 



= J^d^yN{y){detBiy))y^ 



+(det*(y))V2 



)(-+tr*-Hy)) 



= / d3yAr(y)(det5(2/))V2[^(l/^r^(y))(det*(y))V2 



-tv^-\y) 



A 



tr* 



-(det*(y))V2(i/*^(y))2l^(3)(x,y). (28) 



The quantities {l/^f{x)) are actually the reciprocals of the respective eigen- 
values ^f, and therefore are defined only for 7^ and likewise for the 
reciprocals of the squares. Upon integration of the delta function, we have 



1^ = iV(detS)V2 [-Vdii*(^,)-^ + li^frVd^{^ + tr*-i) 



(29) 



where we have used the Liebniz rule in conjuction with the nondegeneracy 
of the CD J matrix ^ae- 

Moving on to the variational derivatives with respect to the configuration 
variable we have 



Assuming det-B 7^ wc can include a factor of 1 in the form (dcti?)^^/2(det-B)^/2, 
which allows us to extract the Hamiltonian constraint from the last factor 
of (30). The result is 



6Xf{x) is ^ SXf{x) 

= ^l^dMB-\y))t^^^N{y)H{y)5(^\^,y). (31) 
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The variational derivative with respect to the configuration variables on 
configurations of nondegenerate Ashtekar magnetic fields reduces to 

upon integration of the delta function in (31). Immediately enter quantities 
i?* which have not been defined on the starting phase space [X^ j'^f). In 
order to explicitly evalute (32) we would need to find dB^/dX^ , which 
requires one to be able to express X^ explicitly in terms of the Ashtekar 
magnetic field B^.^ Whether it is possible to do this or not in the full theory 
for our purposes is immaterial, since on configurations of nondegenerate i?* , 
(32) is directly proportional to the Hamiltonian constraint which we will 
require to vanish.^ 

From (29) and (32) we can write the Hamilton's equations of motion on 
the gauge-invariant, diffeomorphism-invariant subspace of general relativity. 
These are given by 

= ^ = ^(^(*/)-'^ - (deti?)V2(det*)V2(^,)-^) = (33) 



and 



which can be seen as having arisen from the first order starting action 

1 = J d^x ^/X-^ - ATVdet^Vdct^l^^ +tr*-i) . (35) 

The evolution of the system must be confined to the constraint surface H = 
as a requirement of consistency of the system. Therefore, all terms pro- 
portional to the Hamiltonian constraint must vanish and the equations of 
motion reduce to 

X/ = -iV(det5)i/2(^g^^)i/2^^^^-2. ^^ = 0. (36) 



^Since the transformation from Ashtekar variables into the new variables (A",aa) — > 
{X-^,'^f) is noncanonical, then it may be nontrivial to do this in closed form in the full 
theory. 

^As part of the requirements of a consistent theory of constrained systems, the con- 
straints surface must be preserved by Hamiltonian evolution and vice versa. 
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Note that when the Hamiltonian constraint H is satisfied, that = 0, 
which means that ^ f can at most be an arbitrary function of spatial position 
X G S. Hence we have that = Xf are independent of time, and the 
equations of motion for Xf reduce to 

Xf = -iV(detS)V2(detA)i/2(A/)-2. (37) 

The time evolution of X^ is completely determined by the time evolution 
of detB, which naively appears to be unspecified based on the starting ac- 
tion (35) since the Ashtekar magnetic field is not explictly part of the 
canonical structure (25). 



3.1 Algebra of constraints 

Some additional variables of this type will include various quantities related 
to the Ashtekar connection 



SX"^ = BiSAf; dX''3 = OfdXf. (38) 

Suffice it for now to use these definitions (38) to make the result of (30) 
more readable. Continuing with (30), we have 



( SX^■■yiz) ^ ( SB:{y) ^ . . 

= \iyi{{B-')1{B-yp'j>NH). (39) 

The derivative in D^^j^ acts on all terms in brackets, and arises from an 
integration by parts. 

We now compute the Poisson bracket of two Hamiltonian constraints 
using 

{H[NIH[M]} - jj ^(^^^^^ - ^^^^J 

{^D%{{B-')-{B-yp''fMH) -M^n] (40) 
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(39). As one normally expects from the Hamiltonian constraint, which is 
a scalar, only the terms with derivatives acting on the lapses M and N 
survive, yielding 

{H[N],H[M]} = U'^iNdkM - MdkN], (41) 
where we have defined 11^ by 

Hfc = y'^''{B-^)^{B-y.OY(^^{^f )-^H - (det5)i/2(det*)i/2(^^)-2^i/.(42) 

The Poisson of two Hamiltonian constraints is proportional to a Hamiltonian 
constraint with structure functions, which implies a closure of the algebra 
and therefore Dirac consistency. Therefore on the space of solutions to the 
Gauss' law and diffeomoeprhism constraints, the dynamics of the Hamilto- 
nian is self-consistent.^ The conclusion is that in the physical Soo variables, 
the algebra of the Hamiltonian constraint by itself forms a first class con- 
strained system [?].^ If one were to rather attempt to interpret (41) as a 
diffeomorphism, then it would imply a nonclosure of the algebra. 

3.2 A toy model 

It is obvious that on the configruation space G F, the magnetic field 
i?^ contains seven unphysical degrees of freedom. While this may be the 
case, the possibility exists that the physical degrees of freedom inherent 
in which are of relevance to the new variables may reside entirely 
within dcti?.^ This notion can be motivated by consideration of the spatially 
homogeneous sector of configuration space F as a toy model. In this case 
the Ashtekar magnetic field is given by 

Bi = e'i%Al + \e'i'^UcA)Al = {deiA){A-^t. (43) 

''The algebra is still an 'open algebra' in the sense that there exist structure functions 
dependent on the phase space variables. 

®In the original Ashtekar variables, the algebra of the Hamiltonian constraint yields a 
diffeomorphism constant and therefore does not close. Prom the vantage point of the Soo 
variables, the Hamiltonian constraint would form its own subalgebra, with the diffeomor- 
phisms constituting a set of second class constraints. 

®The relation {X^,X'^,X^) ~ (61,62,63), where ha are the eigenvalues of the Ashtekar 
magnetic field seems appealing since this would lead to the relation detS = 616263. A 
rigorous proof of this is beyond the scope of the present paper, therefore we will for the 
time being regard it as a conjecture. 

^°We must emphasize that the homogeneous sector of F does not necessarily mean that 
we are considering minisuperspace. Even if the configuration space variables are spatially 
homogeneous, this doncs not necessarily imply that the spatial homogeneity automatically 
extends to the momentum space variables 
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In (43) we have used the fact that the structure constants fabc = ^abc s^re 
numerically the same as the epsilon symbol €ijk, a property unique to three 
dimensional space S. While X-^ may not have a simple expression in terms 
of the Ashtekar variables, the quantity T = + + X^ integrated over 
S yields 



J SxT{x) = J J Bi{x)5A'}ix) = J ti(^AAdA + ^AAAAA^ = Ics[A] (44) 

which is the Chern-Simons functional of the Ashtekar connection Af. In 
the spatially homogeneous sector of configuration space, we have from (43) 
that deti? = (detA)^, and also from the integrand of (44) that T = detA. 
Summing (37) over the index / we have 



3 

t = Y,X^ = -NvxT, (45) 
/=i 

where we have defined 




A1A2 ^ Ai ^ A2 



4AiAi + Ai + A2 



(46) 



using (20). One can immediately deteminc the evolution of deti? in time, 
which fixes the time evolution of the combination T = X^ + X"^ + X^ , by 
integrating (45) to yield 



r(t) =roe-^^^'^(*) (47) 

Note from (47) that there is still freedom in the choice of the lapse function 
N , which determines the manner of evolution normal to spatial hypersur- 
faces S. For numerically constant N this leads to exponentially inflating 
solutions. 



4 Quantization of the unconstrained phase space 

Upon quantization we must first promote the dynamical variables X^ and 
f to quantum operators X^ and ^/ and the Poisson brackets (26) to 
equal-time commutation relations 
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[X^(x,t),*,(y,t)] = G5f6^^\^,y). (48) 

In the functional Schrodinger representation, the quantum wavefunction is a 
functional of the configuration variables ip = and the basic operators 

act respectively by multiplication and by functional differentiation 



Xf{^,t)iP[X] = {Xf{^,t))iP[X]; 
^f{^,mX] = {hG)^j^^^^[X]. (49) 



4.1 Normalizability of the quantum wavefunctional 

There are two main possibilities for Hilbert spaces in the Soo-CDJ variables, 
depending on the signature of spacetime. For spacetimes of Euclidean sig- 
nature, all variables can be regarded as real. Hence the measure in this case 
can be given by 

3 

DfiEucliX) =DX=l[(^ i^-^SXf{yi) (50) 
/=i X 

Upon the rescaling \f ^ iXf the wavefunctions, which are now a pure phase, 
become delta-functional normalizable 

-i{hGy^ [ (fxCf{x)Xf{x) 

= \{K^}{^)-Qi^))- (51) 

X 

The interpretation is that spacetimes of Euclidean signature not globally 
having the same algebraic classification are orthogonal. 

For spacetimes of Lorentzian signature, it is not possible to impose or- 
thogonality by delta functions using the Lebesgue measure (50). Since we 
require square integrable wavefunctions as a necessary condition for a sen- 
sible quantum theory, we must use a measure which falls off faster than 
the wavefunctions can blow up as X — > oo. Hence, we define the physi- 
cal Hilbert space for Lorentzian signature spacetimes as the set of square 
integrable functions ^'Ai-'^] £ Hpi^yg with respect to the Gaussian measure 



exp 



exp 
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DllLor{X) 



X,/ 



-u-^ J (fxX\x)Xf{x) 



(52) 



The inner product of two physical states for Lorentzian signature is now 
given by 



{^x\^dLor = T{ / ^^(°)'^^^exp[-z.-i / d^xX\x)Xf{ 



exp 



{hG)-^ j d^x\){x)x\x) exp {hG)'^ j d^xQ{x)Xf{x) 



{x) 



exp 

exp 



v{hG)-^ / d^x\){x)Q{x) 



(53) 



The requirement of normaHzabihty is equivalent to square integrability of 
\f{x) and Cf{x). On a final note, we point out that the pure Kodama state 
^Kod exists as a member of the physical Hilbert space ip^, namely for the 
choice Ai = A2 = A3 = 
of Petrov type O. 



■|. In this case, the state would label a spacetime 



5 Quantum Hamiltonian constraint: revisited 

There are two main considerations when quantizing the Hamiltonian con- 
straint, namely operator-ordering and regularization. Operator ordering is 
not an issue in our case since all momenta occur to the right upon quantiza- 
tion, which leaves regularization of products of momenta. Let us define the 
quantized version of the Hamiltonian constraint as the result of averaging 
the constraint with all of its transposes, so that it is Hermitian. Hence, upon 
quantization we have that 



= {^1^2 + ^'2^1 + ^2^3 + ^'3*2 + *3^1 + ^1^3) 
-|-A(*il'2*3 + 5'2^'3^1 + ^'3^'l*2 + ^'2^'l^'3 + ^1^'3^'2 + ^3^'2^l) , (54) 

where Ifg and Ifgh are tensors which impose the desired symmetry averaging 
upon contraction with the appropriate factors of ^f. 

To regularize the Hamiltonian constraint wc make use of the results 
of [?], which uses a point-splitting regularization procedure combined with 
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smearing of each individual factor in the operator product. We make use of a 
regulating function /e(x,y), which is a continuous function of a regularizing 
parameter e, such that 



(55) 



for all smooth test functions ip G C°°(E). The ingredients of the quantized 
Hamiltonian constraint are given by 



H„Jx) 



(56) 



In order to obtain a well-defined action of the constraint of the wavefunc- 
tional, it is not necessary to regularize the first operator since it already 
produces a finite result upon action on the wavefunction. To regularize we 
start with the quadratic term, going into the Schrodinger representation 



^^.^H^G)7^ci3,A(x,y)^^,^|^ 



(57) 



Using the general form of the wavefunction as a third-quantized Wilson line 



ip = exp 



J d^x J Xf[X{x)]SXf{ 



x) 



(58) 



the first functional derivative brings down a factor of A j as in 



i^G)j^MX]=\f[X{xmx]. 



(59) 



As previously mentioned, it is not necessary to regularize this first factor 
(59). Continuing from (57) and regularizing the second factor, we have 



{hG) ^ d3y/,(x, y)^^i-^(Ae(x)V') 



HG 



<5(3)(x,y) + Ae(x)Aa(y) 



= ({nGm){^)^+ I dV.(x,y)A„(x)Ae(2/))v. 



dX'^Jx 



(60) 



where we have defined /e(0) = /^(x, x) and we have used (59). 
Moving on to the cubic term, we have 
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H}ae^ = ^dVe(x,y))^dV.(x,y)^/(^)^a(y)^e(x)V 

= {hGf [ d^zM^,y) [ dV6(x,y) 



6Xf{z)dX''{y)5X^{x) 



tj,. (61) 



The first two functional derivatives in (61) have already been evaluated in 
(60). Hence, continuing from (60), we have 



/ dVe(x,y)A„(x)Ae(j/))V' 



= ^dV.(x,z) {hGffM{ 



32 A, 



.dXfdX^ 



+{hG) 1^ dV.(x,y)(A„(:r)(|^) ^5(3) (z, y) + Ae(y) (1^)^5(3) (x,z) 

+ihG)MO)Xfix) (1^) ^ + Xf{z) ^ d3y/,(x, y)A„(x)Ae(y)j ^. (62) 

We now evalute each individual term of (62). The term of order (hG)"^ is 
given by 

The term of order HG has three contributions. The first contribution is given 

by 

{hG [J^ dhM^, z) ^ /,(x, y)Xa{x) (^) ^,5(3) y)" 
= (nG)Aa(x)^dVe(x,y)/.(x,y)(|^)^ ^ (^G)<,^(x) (64) 

where we have integrated d^z to eliminate the delta function. The second 
contribution of order KG is given by 



(hG) [J^ d^a^, z) ^ d'vM^, y)Ae(y) (|^)/^'^ (x, z) 
(nG/,(0))(|^)^^dV.(x,y)Ae(y) ^ (nG/,(0))/3^,;(a:) 



(65) 
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where we have integrated d^z to get rid of the delta function, in conjvmction 
with factorization of the x dependence. The third contribution of order hG 
is given by 



= hGfM 



' ^ ' dhM^,z)\fiz) = [hGfMHaefi^)- 



(66) 



The semiclassical term is given by 



s: 



aef 



X) 



Xa{x) d^zU^, z)Xjiz)) d'zU^, y)Ae(y) 



(67) 



The next step of the regularization procedure is to isolate the singularities in 
the constraint to poles in e, in conjunction with taking the e ^ limit. To 
this end we will make use of a technique widely used in field theory, namely 
that it is safe to replace quantities with their limits as long as they do not 
blow up. We will repeatedly make use of (55) in the e — limit for all 
quantities occuring under an integral sign, isolating any poles to the orders 
of singularity. This is necessary only for the terms of order hG. 



.o<e/(^) = A„(x)lim,^o^ A/e(x,y)(/e(x,y)(|^) j = /,(0)A„(x)(|^)^(68) 

where we have identified the term in (68) in brackets with ip G C°°{T,). 
Next, we have 

lim,^of3aef{x) = (^^)^lime^o^c?3y/e(x,y)A3(y) = (^^)^^e(^) (69) 
where we have identified <^(y) = \e{y)- Next, we have 



lim^ 



^07ae/(^) = (|^)^lim,^o l^d'yfe{^,z)Xfiz) = Xf{x) (70) 

where we have identified (p{z) = Xf{z). For the semiclassical term we have 



lime-^oSl^f{x) = Xf{x)Xa{x)Xe{x) 



(71) 



The coefficients of the factors of HG and (hG)'^ are finite, which enables the 
isolation of the poles to powers of /e(0). 
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Collecting all of the contributions, we contract the ingredients of the 
quantized Hamiltonian constraint with symmetry tensors lae and Ijae which 
impose the required symmetry (54), obtaining a regularized Hamiltonian 
constraint of 



+ 



1^.(72) 



The action of the regularized Hamiltonian constraint in the small e limit can 
be written in the form 



H'{x)iP = {qo{x) + {hGfM)<li{x) + {hGfMfQ2{x))'il^. (73) 



5.1 Physical Hilbert space 

In the usual methods of loop quantum gravity, one chooses a kinematical 
Hilbert space iJfcjn at the level prior to implementation of the constraints. 
In the Soo/CDJ variables one may choose wavefunctionals of the form 



= exp 



(74) 



where A/[X] are at this stage arbitrary functions of X. Equation (74) al- 
ready satisfies the diffeomorphism and the Gauss' law constraints, which 
leaves remaining the quantized Hamiltonian constraint. Therefore, we will 
need to apply (73) in the limit of removal of the regulator. Hence, 



= iim,^o(go(x) + {hCfMhii^) + {hGMo)fq2{x))'ip. (75) 

Note that as e ^ 0, we have that /,(0) ^ d^^\0) and (/e(0))2 {6^^\0)f. 
In fact, these delta functions of zero are precisely what one would obtain as 
the numerical coefficients of gi and q2 if one were to formally compute the 
action of the Hamiltonian constraint unregularized on Since the Hamil- 
tonian constraint must be satisfied independently at each point, a sufficient 
condition is that qo = qi = q2 = 0. We require that all wavefunctions must 
be independent of the specific form of the regulating function and must 
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be independent of any regulating parameters in the limit when the regulator 
is removed. ^-"^ Since e is arbitrary, it follows that qi = q2 = is necessary 
since these terms must vanish faster than any possible regulator can blow 
up, in order to satisfy the constraint. The result is a set of three equations 
which must be satisfied, namely 



0; 

0; 

0. (76) 

Equation (76) is a set of three equations in three unknowns Ai[X], A2[^] 
and A3[X], and we must construct a Hilbert space of states of the form 
(74) solving these equations. The general solution involves an integration 
in the functional space of fields T, which involves constants of functional 
integration with respect to at each point at which the Hamiltonian 
constraint is satisfied. We will be content, for the purposes of making 
contact with the previous sections, to consider wavefunctions for which Xj 
is independent of the dynamical variable . In this case the second and 
third equations of (76) automatically vanish, and we are left with the first 
equation which is the semiclassical part of the Hamiltonian constraint 

2(AiA2 + A2A3 + A3A1) + AA1A2A3 = Vx. (77) 
The general solution to (77) is given by 



go = /"'A„Ae + AF'^^A/A„A, 



dX"- dXf 



A3 = -(xTT^i^). (78) 



^AiA2 + Ai + A2 



The corresponding wavefunctional, which then acquires the labels Ai and 
A2, is given by 



V'A[^]=exp (hG) J (fx\f{x)Xf{x) 



(79) 



Since the configuration variables X-^ are now conjugate configuration vari- 
ables for the eigenvalues of the symmetric part of the CD J matrix ^ae, then 
there is no dependence upon the variable conjugate to ipd- Hence we have 



^^This is the standard required of a finite state of quantum gravity by our interpretation. 
^^The integration over T therefore introduces continuous functions with each quadrature, 



which can be used to label the states. 
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that i/^x £ Ker{H (~]Hi}, and whatever that conjugate variable is, it must 
be unphysical. Additionally, ^ -f^erjGa} since Ai and A2 fix the required 
iS'0(3, C) frame for a given B^. 



ipx[X] = expf / d^xXf{x)Xf (x) 



(80) 



The functional (80), which satisfies the quantum Hamiltonian constraint, 
also satisfies the diffeomorphism constraint. Since the configuration vari- 
ables are now conjugate configuration variables for the eigenvalues of 
the symmetric part of the CD J matrix ^ae^ then there is no dependence 
upon the variable conjugate to ijjd- Hence we have that ij^^^ G Ker{H f]Hi}, 
and whatever that conjugate variable is, it must be unphysical. 

We will ultimately move on to the Gauss' law constraint, but let us first 
address the normalizability of the wavefunctions (80) . Since in the worst case 
the variables X^ may be complex, we require Hxin C L^(r, D/XLor(-'^)) to 
be the set of square integrable functionals on the functional manifold F, with 
Gaussian measure 



r>AiLor(^) =-DXe-^^(^)^^(^) = JJ(5X(x)exp - j d^xX f{x)xf {x) .(81) 



6 Discussion: Hamilton— Jacobi functional 

Prom (20) one may construct the Hamilton-Jacobi functional 



S= I d^x 



X^X' + MX^ - ( ^ )X 



(82) 



which exponentiated in units of KG yields a wavefunctional 



Va[^] = e^^^ (83) 

In (83) A/ = (Ai , A2 , A3) are a set of labels and the variables X-f = {X^ ,X'^,X^) 
are a set of configuration space variables whose role will become clear. 

To put (82) into a more physically intuitive form, let us rewrite it in 
terms of the trace and the traceless parts. First make the substitution 
X^ = T-X^ -X^ into (82) to eliminate X^ in favor of the trace T. The 
integrand of (82) reduces to the form aX + (3Y — ET, where we have defined 
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X = and Y = X'^. The following nonlinear relation ensues amongst a, 
/? and Ai and A2 



^AiA2 + Ai + A2 



2A1A2 



-9A1A2 + Ai + A2 



(84) 



and 



/3 — a = A2 — Ai — ^ A2 = Ai + /? = a. 



(85) 



Substitution of (85) into (84) yields the following cubic polynomial equation 
for Ai 



(|ai + l)(Ai +(i-af + (2Ai - /3(^Ai + 1)) (Ai + ^ - a) - = 0, (86) 



which simplifies to 



(Ai)3 + Q + /3 - 2«) (Ai)2 - (x + - ")) ^1 - a"^^ - «) = (87) 

To solve (87) using trigonometric functions, it helps to put it into dimen- 
sionless form using the rescaling Ai = ^ , where 7 is a numerical constant of 
dimension [7] = 2 and p is dimensionless. Then (86) reduces to 



/ + V + + C = 0, 



(88) 



where we have defined 



^ = 7'(| + /3-2a); B = + a{P - a)); C = -j^[jaiP - a)) .(89) 
The general solution for the first eigenvalue Ai is given by 



Ai = Ai(a,/3) = -^ + y^^(s-^)r[/3 



-4(C-^ + ^) 



3/2 



(90) 



where r = ^ labels the the three roots of the cubic polynomial for n 
(0, ±1), and we have defined the Chebyshev polynomial 



T^{x) = sin[asin + r] . 



(91) 
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The solution for the second eigenvalue A2 automatically follows 



A2 = A2(a,/3) = Ai(a,/?)+/3-a. (92) 
Making the definition 

4 Ai (a, /3)A2 (a, /3) + Ai (a, /3) + A2(a, /3) 
one can construct a wavefunction by exponentiating (82) of the form 

^l}^p[X,Y,T]=e^^\-{hG)-^ j (i^x(a{x)X{x) + (3{x)Y{x) - E^p{x)T{x)) 

Equation (94) is a gauge-invariant, diffeomorphism invariant wavefunctional 
which satsifies the classical Hamiltonian constraint by construction. As an 
aside, for a = /3 = 0, we have that £^00 = £ (^^) ^'educes to the pure 
Kodama state 

^Kod = ^oo[X,Y,T] = e-6(^A)-^/.^. (95) 
Note that the trace T = + + is directly related via 

T = J I B'JA^ = j ti(^A AdA+'^AAAAA)j, (96) 

which is the Chern-Simons functional for the Ashtckar connection. The 
exponential of (96) yields the Kodama state ^ Kod-, which is the both a 
semiclassical manifestation of DeSitter spacetime as well as a quantum state. 
For more general case (94), the semiclassical orbits encompass more general 
solutions to the Einstein's equations. 
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